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ABSTRACT 

We find significant fluctuations of angular momentum within the convective helium shell of a pre-collapse 
massive star - a core-collapse supernova progenitor - which may facilitate the formation of accretion disks and 
jets that can explode the star. The convective flow in our model of an evolved Mzams = ISMq star, computed 
with the sub-sonic hydrodynamic solver MAESTRO, contains entire shells with net angular momentum in 
different directions. This phenomenon may have important implications for the late evolutionary stages of 
massive stars, and for the dynamics of core-collapse. 
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1. INTRODUCTION 

Theory suggests that massive enough stars develop a core 
of iron-group elements, which upon reaching a critical mass 
undergoes a catastrophic collapse followed by an energetic 
and luminous explosion termed a core-collapse supernova 
(CCSN). Although stars are close to being spherical objects, 
there are phenomena for which a spherically symmetric one¬ 
dimensional treatment may be insufficient, both in the nearly 
hydrostatic states during the long stellar lifetime prior to col¬ 
lapse, and in the hydrodynamic explosion itself. 

In stellar evolution modeling, convective mixing and trans¬ 
fer o f energy is a nota bly multi-dimensional phenomenon 
(e.g., iBeeck et'HIUOlSl) . Convection is usually treated us¬ 
ing the approximate semi-empirical Mixing-Length Theory 
(MLT). This simplified treatment might be inadequate in some 
cases. Ongoing studies aim at improving the theoretical mod- 
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it might be important to look at multi-dimensional 
effects such as the standing accretion-shock insta¬ 
bility (SASI; e.g., [ Blondin. Mezza cappa & DeM arinol 
2003t i Blo ndin & Mezzacappal 1200^ iFernandezi 1201 dk 
Papish, Nordhaus & Sokeil (120151) ^ggest that neutrino- 
driven mechanisms might not work. 

The multi-dimensional phenomena in evolution and ex- 
plosion modeling may be interconnected. Rec ent studies 
(ICouch & Ottl|20lll20lHMueller & Jan^lMsh found that 
non-spherical perturbations in the progenitor may increase 
post-collapse turbulent pressure, significantly influencing the 
dynamic behavior in the shocked region surrounding the 
forming NS. In these studies, large velocity fluctuations were 
introduced in the pre -collapse core. In an earlier paper 
dGilkis & SokerllTolSh we found that such velocity perturba¬ 
tions imply sufficiently large specific angular momentum fluc¬ 
tuations for the formation of intermittent accretion disks and 
a jittering-jets explosion. 

Among the alt ernatives to the neutrino-driven explosion 
mechanism (see iJankal 120121 for a review), most promi¬ 


nent are jet-driven explosions (e.g.,[ .^eBlanc & Wilsonlll970t 
iKhokhlov et al.1 Il999t iLazzati et ak 1201 21 . Another inter¬ 
esting m echanism is the collap s e-induced therm onuclear ex¬ 
plosion (iKushnir & Katzl 120151: lKushnirll2015allbh . The jet- 
driven scenario is inherently multi-dimensional, where a 
pre-collapse rapid core rotation may facilitate t he formation 
of ma gnetorotationally driven bipolar outflows (iMdsta et al.l 
l2015h . We note that the above cited works require a rapidly 
rotating pre-collapse core, and hence limited to a small frac¬ 
tion of all CCSNe. 

A different approach of jet-driven explosions tries to ex¬ 
plain all energe tic CCSNe, for cases of rapid as we l l as slo w 
stellar rotation (ISokeriUOlft iPapish & Sokedl201 IL l2012allR 
l2014allbll . In this model accretion from material with varying 
angular momentum in a collapsing star leads to the formation 
of accretion disks with varying axis direction. These inter¬ 
mittent accretion disks launch jets with different directions, 
termed ‘jittering jets’, that might explode the star with the ob¬ 
served energy. The jittering jets model faces the challenge 
of supplying the required angular momentum for an accre¬ 
tion disk, and consequently forming jets. The required an¬ 
gular momentum for jet formation may be lower than usu¬ 
ally implied, when a d ynamo operates in a shear layer around 
the pro to-neutron star (lAkivama et al.l200^ISchreier & Sok^ 


Recent studies have concentrated on non-spherical pertur¬ 
bations in the close vicinity of the iron core just prior to 
collapse (e.g., ICouch et al.l 120151) . We look a bit earlier in 
the stellar evolution and further out in the star for a possi¬ 
ble source of angular momentum. While the outer hydrogen 
region probably has tremendous fluctuations of angular mo¬ 
mentum, the sudden change in the gravitational potential as 
a result of neutrino emission d uring NS formation, results i n 
the expulsion of the envelope (iLovegrove & Wooslevll20l'^ . 
Analytic estimations show that the convective helium shell of 
an evolved massive star may have large fluctu ations of angu¬ 
lar momentum as well (iGilkis & Sokeil 120141) . In a scenario 
where after core-collapse the inner silicon/oxygen regions fail 
to produce a sufficient driving force for a supernova explo¬ 
sion, accretion from the helium region might result in the 
formation of intermittent accretion disks and jets. Accretion 
from material with varying angular momentum in a collaps¬ 
ing star can lead to an energetic supernova explosion in the 
jittering jets model. 

In the present study we explore the properties of the pre- 
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Fig. 1.— Detailed composition for our stellar model of a non-rotating 
JW^ZAMS = 15M0 star with metallicity of Z = 0.014. The profile is 
presented at an age of f = 12.5 Myr, about 5.5 yr before explosion. The 
stellar parameters at this stage are as follows. Stellar mass of M = IS.SATq, 
effective temperature of T = 3235Ff, photospheric radius of R = 902i?Q 
and luminosity of L = 8 X IO^Lq. The convective helium shell simulated 
here extends from rin = 3.7 X 10® cm to rout = 2.8 X 10^® cm. Results 
below will all be presented at this time and for this model. 

collapse convective helium region. We perform a three- 
dimensional hydrodynamic simulation of the convective flow, 
strengthening our preliminary analytic estimations. The he¬ 
lium region does not change much in the relatively short time 
left until collapse, so that our results are relevant for the final 
stage of the stellar evolution. In section|2]we describe the nu¬ 
merical setup and method. Our main results are presented in 
section [3 Some implications for core-collapse are concisely 
discussed in section0] We summarize in section|5] 

2. NUMERICAL SETUP 

We use a stellar model constructed by Modules for Ex- 
periments in Stellar A strophysics (MESA version 5819; 
iPaxton et al.ll201 ll 120131) . with an initial mass of Mzams = 
15M0 and metallicity of Z = 0.014. Due to stellar 
winds calculated here with the so-called ‘Dutch’ schern e (e.g., 
iNugis & Lamersl2000l:iVink, de Koter & Lamersl2001h the fi¬ 
nal mass is IS.SMq. Effects of rotation and magnetic fields 
are neglected. We examine the star during the helium shell 
burning stage, about 5.5 years prior to collapse. At this 
stage, the star is a red supergiant (RSG) with an effective 
surface temperature of T = 3235iT, photospheric radius of 
R = 9O2i?0 and luminosity of L = 8 x 1O‘*L0. The detailed 
composition of the stellar model is shown in Eigure[T] 

There are two reasons for simulating the star at this time. 
(i) The convective helium shell may be of interest as a po¬ 
tential supply of material with angular momentum for accre¬ 
tion onto the newly-formed compact object following core¬ 
collapse, in a scenario where the inner shells fail in driving 
a successful explosion. The convective hydrogen envelope 
is less relevant, as it is likely to be expelled due t o decrease 
of the NS gravitational mass by ne utrino cooling (iNadezhinI 
119801: iLovegrove & WooslevI l2013h . The structure outside 
ro-core ~ 3.5 X 10® Cm remains essentially unchanged un¬ 
til core-collapse, so the results presented here are relevant for 
the core-collapse epoch. If we had taken an earlier point of the 
stellar evolution, the different core structure would have less¬ 
ened the pertinence of our simulation, (ii) While in the final 


stage before core-collapse there are several convective shells 
(silicon and oxygen, in addition to helium and hydrogen at 
the simulated stage), the complex structure would have made 
the simulation more difficult to run. We hope the methods 
detailed here will inspire future works on angular momentum 
fluctuations in the silicon and oxygen shells in core-collapse 
supernova progenitors. 

Convection in MESA is according to the Mixing-Length 
Theory (MET), which gives for our model a convective outer 
hydrogen envelope, and a convective helium region surround¬ 
ing a radiative oxygen core. The outer part of the helium re¬ 
gion is radiative as well. In this study we are not interested 
in the outer envelope, and focus on the convective helium re¬ 
gion, extending from r-m = 3.7 x 10® cm and M-^ = 2AMq, 
to Tout = 2.8 X 10^® cm and Mout = If and when 

this region collapses to the core, the free fall time from the in¬ 
ner and outer boundaries of the helium convective regions are 
= 14 s, and fff_out = 228 s, respectively. Interestingly, 
this range is compatible with the duration of long gamma ray 
bursts (LGRB). This is not the case studied here, but the be¬ 
havior of the helium convective layer might be relevant to the 
formation of jets in GRBs. 

The inner region of the one-dimensional MESA model is 
mapped into a three-dimensional grid in the lo w Mach num¬ 
ber s tellar hydrodynamics code MAESTRO (iNonaka et al.l 
I 2 OIOI) . The MAESTRO code was developed mainly for the 
purpose of studying convecti on in white dwarfs exp eri encing 
thermonuclear defl agration (lAlmgren et al.l l2006allbl 120081 
IZingale et al.l 1200^ . though more recently als o for model¬ 
ing convection preceding Type I X-ray Bursts (iMalone et al.l 
I 2 OI IL 120141 IZingale et al.l 120151) and conyection in massive 
stars (iGilet et al.l 2013h . This latter study motiyated us to use 
MAESTRO for our purpose - the study of angular momen¬ 
tum fluctuations in a conyectiye region of a massiye star. Eor 
this purpose we need a full sphere domain, with no imposed 
reflectiye boun daries for reductio n of computational costs. 

Similarly to iGilet et akl (120131) . we initialize the yelocity 
held with small perturbations (much lower in magnitude than 
the expected, and obtained, flow yelocities), and use yelocity 
damping near the edges of the simulation domain (far from 
the regions of interest). We further added yelocity damping 
in the oxygen core, as we encountered numerical problems 
in this region of the simulation - spurious yelocities and heat¬ 
ing. This may be due to the absence of radiatiye transfer in the 
simulation (important in the non-conyectiye core) or because 
of insufficient resolution. We ran simulations with damping of 
the entire core, part of it, and no damping at all. The simulated 
conyectiye helium region changed little between the simula¬ 
tions, and the inner oxygen core is of no particular interest in 
the present study. We therefore conclude that our results are 
robust, although we may need to address this issue in a future 
work. 

We haye incorporated an approximate reaction network in 
the code, which includes only the triple-alpha process (the 
only releyant process in the re gion of interest). The energy 
generation rate is according to iKippenhahn, Weigert & WeisI 
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e3a = 6.272p2x| • (1 + 0.00158rg-° ®®) 

X [2.43 X exp (^-13.490rg'^/^ - (Tg/O.lS)^) 

X (1 + 74.5T9) + 6.09 X lO^Tg”^/^ exp (-I.O54/T9)] 

X [2.76 X 10^Tg"^/^exp (-23.570Tg”^/^ - (r9/0.4)^) 

X (1 + 5.47r9 + 326r|) + 130.7Tg”^/^exp {-3.358/Tq) 

+ 2.51 X 10'‘Tg"^/^exp(-20.307/T9)], (1) 

where £3^ is in erg gr“^ s“^, p is the density, X 4 is the 
helium mass fraction, and T9 is the temperature divided by 
10® K. We have neglected composition changes, as the 
nuclear timescale is significantly longer than the dynamical 
timescale. 

The Adaptive Mesh Refinement (AMR) capabilities of 
MAESTRO were employed to keep a grid resolution better 
than O.lTfp throughout the convective helium region, where 
Hp is the pressure scale height. We performed several sim¬ 
ulations with lower finest resolution (higher resolutions are 
prohibited with our current computational resources). The 
magnitude of the angular momentum, which was our primary 
concern, is similar in the different simulations. The maximal 
Mach number of the convective motion and the average ve¬ 
locity, however, increase with decreasing resolution. Clearly, 
simulations with higher resolution are required to fully under¬ 
stand the phenomena investigated here, but for this introduc¬ 
tory study our results are sufficient. A detailed analysis of the 
lower-resolution full simulations is presented in Appendix lAl 
and additional simulations of one eighth of the domain and 
higher resolution are presented in AppendixlBl 

3. SIMULATION RESULTS 
3.1. Angular momentum 

We ran the simulation for a physical time of close to f = 
20000 s. This is about 85 times the free fall time from the 
outer convective boundary fff_out, and about 8 times the typ¬ 
ical convective overturn time fconv = Dconv/vconv ~ 2500 s. 
The latter is calculated from the the typical diameter of the 
vortices Dconv ~ 5 x 10® cm, and the typical velocity in the 
vortices of Uconv ~ 2 x 10® cm s“^, as seen in figures|2]and 
[3 presented at the end of the simulation. Eigure |2| shows the 
flow morphology, with color indicating the local specific an¬ 
gular momentum calculated from the flow velocity. Eigure [3] 
shows a close-up of a region from Eigure [2] with color indi¬ 
cating the magnitude and sense of the velocity, to emphasize 
the large-scale vortices that develop in the core. 

Eigure m shows the evolution with time of the global com¬ 
putational diagnostics of the maximal Mach number and the 
averaged velocity magnitude. It can be seen that after leaving 
the initial close-to-zero velocity setup, a steady flow devel¬ 
ops after a time of about fconv Another important conclusion 
drawn from this figure is that the typical convective velocities 
are much higher than values given by the mixing-length the¬ 
ory. Similar results o f higher than the M ET velocities have 
been obtained by, e.g.. lArnett et. all(l2009h . This result holds 
even when a velocity lower by a factor of 2 is taken, as might 
be the case with higher resolution (see Appendix iBli. 

Eigure [5] shows the evolution of the angular momentum in 
several shells. We notice the significant net angular momen- 



Fig. 2.— Specific angular momentum in the xy-plane at t = 20000 s 
in the simulation. Blue-shaded regions have a clockwise azimuthal velocity 
component (jz < 0), while red-shaded regions have a counterclockwise az¬ 
imuthal velocity component (jz > 0). White regions have no azimuthal ve¬ 
locity component (jz = 0). The color-coding runs from —2 x 10^® cm^ 
up to 2 X 10^® cm^ s~^, although the specific angular momentum exceeds 
this value at some points. The inner oxygen core is omitted from the figure. 



Fig . 3.— Azimuthal velocity color-map for a close-up of a region from Fig¬ 
ure with added velocity arrows for illustration of the vortices in the flow. 
The color-coding runs from an azimuthal velocity of —2 x 10^ cm 
up to 2 X 10® cm s“^. Similarly to Figure[2] blue-shaded and red-shaded 
regions have clockwise and counterclockwise velocity components, respec¬ 
tively. The velocity arrows are scaled - larger arrows represent higher ve¬ 
locities (an arrow spanning 10® cm on the plot, for example, represents a 
velocity of 2 x 10® cm s~^). 
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Fig. 4.— Left: Maximal Mach number in the simulation domain. The val¬ 
ues reached in the simulation, of M ~ 0.15, are two orders of magnitude 
above the mixing-length theory value for this model (AImlt ~ 10“®). 
Right: Average velocity in the simulation domain. This too is much higher 
than the velocities of the MLT, which are around umlt ~ 10® — 2 X 
10® cm s“® for our stellar model. The disparity to the MLT is smaller for the 
averaged velocity, as there is an extensive range of velocities in the flow field, 
and the maximal velocities (and Mach numbers) reach values significantly 
higher than the mean. 


turn that is developed in the shells, and the varying direction of 
the angular momentum axis. In Figure|6]we show the specific 
angular momentum profile of the material within the inner 
part of the convective region at f = 20000 s in the simulation. 
In plotting Figure|6]we used shells of width Ar = 2000 km, 
much thinner than the Ar = 40000 km shells used in draw¬ 
ing Figure |5] The difference in free fall times between the 
two boundaries of the thin shells isf// « 1 — 3 s, equiv¬ 
alent to thousands of dynamical times on the surface of the 
newly formed NS (or BH). Shorter time scales may be of in¬ 
terest for the post-collapse dynamics, but a higher resolution 
is required, accordingly. 

From Figure |4] we note that the average convective velocity 
increases rapidly within the first 2 x 10^ s. The value then 
decreases and increases back, followed by a decrease in the 
time interval 8 x 10^ — 1.8 x 10^ s, after which it increases 
until we terminate the simulation. Variations over a timescale 
ofRiO.5 — IxlO'^s are seen also in the values of the angular 
momentum in thick shells as presented in Figure |5] Both fig¬ 
ures |4] and |5] seem to suggest that the simulation has reached 
a more or less steady state, but with significant fluctuations. 
It would be useful to extend the calculations to check that 
this is actually the case, although this is computationally ex¬ 
pensive. The large amplitudes and the long timescale of the 
fluctuations are the source of the relatively large value of the 
stochastic specific angular momentum of the accreted gas. 

Another indication that we have reached a more or less 
steady state is presented in Figure |6] We see that the spe¬ 
cific angular momentum components in thin shells reach max¬ 
imum values that are about equal to those of the standard de¬ 
viation of specific a ngular momentum calcu lated analytically 
from the MLT as in iGilkis & Sokerl (l2Q14l) . This estimation 
depends on th e shell width, Ar, whi ch in the present work is 
larger than in IGilkis & Soked (120141) . In our previous work, 
the thin shells had a width equivalent to a difference of free- 
fall times that is a small multiple of the Keplerian time for 
a NS, which explains the higher values obtained for the stan¬ 
dard deviation of specific angular momentum there. Addition¬ 


ally, the analytic estimation using MLT values for the present 
work is somewhat problematic, as the flow velocity is higher 
in the numerical simulation than in MLT. On the other hand, 
the length scale taken may need to be smaller than the mixing- 
length, which is close to the scale of vortices here, while the 
relevant length might be only part of a vortex. These fac¬ 
tors compensate each other to some extent, so that the ana¬ 
lytic estimation is close to the fluctuations of specific angu¬ 
lar momentum in the simulation. This argument is an order- 
of-magnitude one, but nonetheless supports the claim that a 
steady state has been achieved. The diameter of the vortex, 
for example, is taken to be of the order of the mixing-length 
(or the radius of the vortex Oc about equal to half the mixing- 
length). This order of magnitude estimate, as presented in 
Figure|6] further demonstrates the merit of the analytic model. 


3.2. Convective luminosity 

To further reveal the convective properties of the 3D sim¬ 
ulations we study now the energy transport property of the 
convective motion. We start by writing the three equations as- 
su med by the MLT. I n doin g so we follow the nomenclature 
of iKippenhahn et ^ (l2012h . The convective flux is given by 
the heat capacity per unit mass of the gas Cp multiplied by the 
mass flux pvc, where Vc is the convective velocity according 
to the MLT, and by the temperature difference DT between 
the starting and end point of the convective element 


Fconv = pVcCpDT. (2) 

The temperature difference is given by 

Dr = (v-Ve)^-^r (3) 

Z ilp 

where V = 9 In T/O In P is the temperature gradient (with 
respect to pressure) in the star, while Ve is the temperature 
gradient of the convective element; Im is the mixing length 
and Hp is the pressure scale height. The velocity of the con¬ 
vective element is given by 

t;2=p^(V-Ve)^, (4) 

where S = — {dinp/dlnT)p-, for an ideal gas <5 = 1. 

We now manipulate these three equations as follows. From 
equations (O and (IDi we derive 


DT = 


AT 

gSlm 


which together with equation (|2|i give 


( 5 ) 


Fr., 


= PVc 


ATcp 
gSlm ' 


(6) 


The convective luminosity, equals to the integration of the 
convective flux over a spherical surface, can be written then 
in the form 

Pconv — A(r)pKE5 (^) 

where A (r) = ATcp/g6lm depends only on the radius. In 
the MLT the kinetic luminosity Pke is given by 


Ake-mlt (t) = Anr'^pvl. (8) 


Let us summarize the important assumptions that went into 
the derivation of equation (0. (i) The convective elements 
perform radial motion, (ii) They are destroyed at a radius that 
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Fig. 5.— Angular momentum components of several shells within the convective region throughout the simulation time. The location of the shells is indicated 
in the inset in the right panel. Several shells possess net angular momentum with large values and in different directions. The masses of the shells from the 
inner shell to the outer are 0.25, 0.25, 0.23, 0.22, 0.19 and O.ITMq. The core mass inner to the helium shell is Mjn = 2AMq,. The average specific angular 
momentum is therefore approximately few x 10^^ cm^ s“^, which is lower than the specific angulai' momentum for a Keplerian orbit around a NS or BH. 
However, the partition into shells with widths of Ar = 40000 km is arbitrary. While we see here the large-scale behavior, our simulation might not capture 
some finer details relevant for post-collapse dynamics, as the differences in the free-fall times between the shell boundaries correspond to thousands of orbital 
times around the newly-formed NS or BH. A more focused analysis of the angular momentum distribution within shells at one time is shown in Figure[^ 
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Fig. 6.— Specific angular momentum components in thins shells of 
Ar = 2000 km aX t = 20000 s in the convective zone in the simu¬ 
lation, between r = 40000 km and r = 280000 km. Black dashed 
lines show an analytic approximation of th e standard devia t ion of specific 
angular momentum in the MLT, following IGilkis & Sokeil <20141) so that 

a (j) = (2/27)^/^ (Ar)“^/^ with ac taken as half the mixing- 

length and Vc the MLT convection velocity. The specific angular momen¬ 
tum for a Keplerian orbit around a newly-formed NS or BH is approxi¬ 
mately jKep ~ few X 10^® cm^ s“^; the values for the thin shells shown 
here is lower, reaching jmax (fluctuations) 6.85 x 10^"^ cm^ 
O.OSjKep- Although the specific angular momentum is lower than the Ke¬ 
plerian, it is not negligible, and a belt-like structure mig ht form around the 
NS or BH (as discussed bv IPapish. Gilkis & Sokeiir2015h . The difference in 
free-fall times for Ar = 2000 km is 1 — 3 s, corresponding to hundreds of 
orbital times around the newly-formed compact object, so there is sufficient 
time to form accretion belts and jets upon core-collapse. 


is larger by from where they were formed. (Hi) They trans¬ 
fer all their extra thermal energy where they are destroyed. 

In the full 3D simulations we find that the motion is com¬ 
posed of vortices, rather than radial motions (Fig. O. This 
implies that the hotter gas that is flowing at the top of a vortex 
loop is not destroyed, hence does not transfer to the surround¬ 
ings all of its thermal energy. However, the net flux should be 
about equal to that in the MLT, as the structure of the star is 
the same, hence same luminosity. To compare the convective 
luminosity in our 3D simulation to that of the MLT, we need 
to replace Lke-mlt given in eqaution (l8]l by 

pTT /'‘ZtT 

iKE-3D (r) = / sinOdO d4>pv'^ {v ■ f). (9) 

Jo Jo 

where the density and velocity now depend on the location, 
p (r, 9, (p) and v (r, 9, p), respectively. To take into account 
the numerical grid, in calculating Lke-3D (r) equation is 
replaced by 


iKE-3D {r) = 

nr-\-Ar pTT p2-k 

X dr' sm9d9 / dppv'^ {v ■ f'). (10) 

Jr Jo Jo 

In Figure [Tj we plot the kinetic luminosity of the MLT ac¬ 
cording to equation (|8]l, in blue dots, and that of the full 
3D simulation according to equation (fTOl i. in green squares. 
Although the kinetic luminosity of the 3D simulation is not 
identical to that of the MLT, it has the same bulk behavior. 
Our results at other times show a kinetic luminosity profile 
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Fig. 7.— Comparing the kinetic luminosity of the MLT according to equa¬ 
tion i) , in blue dots with that of the full 3D simulation according to equation 
Go), in green squares. 

that changes with time. For example, the peak near the inner 
boundary changes slightly, such that the average value of the 
3D results is less than twice that of the MLT near the peak of 
the kinetic luminosity. Overall, Figure [T] makes us confident 
in the convective properties revealed by our 3D simulations. 
Most important, the convective region is composed of large 
vortices that may translate after collapse to local fluctuations 
in the specific angular momentum of the mass accreted by the 
newly born NS or BH. 

Despite our finding that the typical convective velocity in 
the 3D simulation is much larger than the typical velocity in 
the MLT, the convective luminosity of the two cases are very 
similar. The reason is as follows. In the MLT the convec¬ 
tive velocity is calculated from the assumption that hot con¬ 
vective cells deposit all their extra energy after they move a 
radial distance equal to the mixing-length. This implies that 
the calculation of the convective luminosity of the MLT in¬ 
volves only upward moving gas. In the full 3D simulations 
convective cells are not destroyed when they reach maximum 
radius in the convective zone. They rather perform motions 
along vortices, and deposit a small fraction of their extra en¬ 
ergy when they reach the maximum radius in their cyclical 
motion. Therefore, the integral performed in equation (|9ll (or 
eq. fTOl) involves a large flux downward. Namely, the factor 
V ■ f in those equations posses both positive (for upward mo¬ 
tion) and negative (for downward motion) values. The very 
large flux upward minus the very large flux downward gives 
a value of the convective luminosity about equal to that in the 
MLT, as expected if the numerical simulation converges to a 
solution that yields the correct luminosity through the convec¬ 
tive zone. 

4. IMPLICATIONS FOR CORE-COLLAPSE 

The flow structure of the helium presented in section [3 is 
probably similar to that just prior to collapse, and we can 
look at possible implications on the post-collapse dynamics. 
The scale of the vortices is such that the difference between 
the free-fall times onto the newly formed NS or black hole 
(BH), between the edges of a vortex, is about tens of seconds. 
The local specific angular momentum as shown in Figure|2]is 
higher than the Keplerian specific angular momentum for or¬ 
bit around a newly-formed NS. The large vortices imply then 
that on a short time scale of « 0.01 s, equal to several dynam¬ 


ical times on the newly born NS or BH, a small number of 
vortices contribute to the accretion, and only a small fraction 
from a vortex is accreted, contributing a small mass with high 
specific angular momentum. The differences in free-fall times 
between shells with differing angular momentum axes, equiv¬ 
alent to many orbital times close to compact object (where we 
assume jets are launched), allow for sufficient time for forma¬ 
tion of disks and jets. The average value obtained for shells 
as described in Eigure|6]is lower, and may lead to the forma¬ 
tion of a thick accretion disk (or an accretion belt; for fur ther 
discussion of this general eff ect see lGilkis & SokeJl2015l and 
iPanish, Gilkis & Sokeill2015T) . 

Considering jittering jets, which we suggest explode the 
star, we can make the following crude energy estimate. We 
take only shells with the high specific angular momentum to 
form accretion belts and launch jets. Say we take only parts 
with j > 6 X 10^^ cm^ according to Eigure|6j the fluc¬ 
tuations in narrower shells are larger even. The mass in these 
parts is Ri O.IMq. Accretion of this mass onto the NS of a 
mass of ~ 2Mq, or later a BH, can be channeled into jets with 
energy « 10% of the rest mass of the shell. This amounts to 
« 2 X 10®^ erg. To account for a CCSNe it is sufficient that 
only about 10 perc ent of the time the accr etion belt/disk exists 
and launches jets (IPanish & Sokeill201 11) . We emphasize that 
this last calculation is approximate and somewhat speculative. 
However, it is worth consideration as our main motivation for 
studying angular momentum fluctuations in the helium shell 
was implications for jet-driven explosion mechanisms. Due 
to the high resolution required, the need to explore convec¬ 
tive regions closer to the iron-core, and the uncertainties of 
the jet-formation mechanism, we cannot yet reach definitive 
conclusions regarding the implications of convective regions 
for the jittering jets model. 


5. DISCUSSION AND SUMMARY 

We have presented a novel study of angular momentum 
fluctuations in stellar convective regions, using the sub-sonic 
solver MAESTRO. Further refinement and development of 
this method may enable the study of such fluctuations in dif¬ 
ferent stellar models at different stages (e.g., just prior to core¬ 
collapse). These fluctuations may have implications on the 
core-collapse explosion mechanism, but may als o be relevant 
for th e process of stochastic spin-up of the core (iFuller et alJ 
l2015h . In the present study we found large-scale fluctua¬ 
tions of angular momentum in the convective helium zone 
of an evolved massive star, while the energy transport is 
similar to that expected from MLT, strengthening the valid¬ 
ity of our simulation. Future studies will explore similar 
phenomena in silic on and oxygen shells (as done by, e.g., 
Mueller e t al.ll2016 ), and effects of rotation (as done by, e.g.. 


Chatzopoulos et al 


IMl). 


In this study we focused on the convective helium zone 
surrounding an oxygen core. During the remainder of 
the stellar evolution time, the late burning stages eventu¬ 
ally produce an iron core surrounded by silicon and oxy¬ 
gen shells. We assume that the structure of the helium 
shell does not change much in the relatively short time 
left to collapse. This implies that our results are relevant 
for the final stage of the star, although there may be late- 
stage processes beyond our one-dimensio nal stellar modeling, 
such as significant wave-driven mass-loss (lOuataert & Shiode 
2012 JShiode & Ouataerll2014t) or expansion (iMclev & Sokei 


12014. 


The mass inner to the convective region studied here is 
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Min ~ 2AMq, and a jet-driven explosion powered by ac¬ 
cretion from the helium zone will leave a heavy NS. Loss 
of mass-energy to escaping neutrinos may reduce the rem¬ 
nant mass to around Mns ~ 2Mq. For a lower mass 
NS to form (e.g., with Mns ~ 1-4M0), a jet-driven ex¬ 
plosion will have to be powered by accretion from the sili¬ 
con/oxygen shells. For more massive stars, the helium shell 
may surround a core massive enough to form a BH upon col¬ 
lapse. If the accretion of helium pow ers an explosion, then 
the relevance of ‘failed supernovae’ (iKochanek et all 120081: 
lO’Connor & Otlll201 ll iKochaneldlm^ or ‘very weak su- 
pernovae’ (iLovegrove & Wooslev 2013h may be lessened. A 
supernova leaving a BH remnant may be prevalent for very 
massive stars - this may b e the case for, e.g., SN 2005gl 
(iGal-Yam & Leon^l2009l) . 

We can consider the present results on the source of 
stochastic accreted angular momentum from the helium 
shell in combination with our previous results on a pos¬ 
sible stochasti c angular momentum accretion from the 
silicon shell (iGilkis & Soked 1201 5h , and the results of 
iPapish, Gilkis & Sokerl (1201 5h regarding the stochastic angu¬ 
lar momentum accretion as a result of the standing accretion- 
shock instability (SASI). We then come at an emerging idea 


where the newly formed NS or BH accretes mass with notable 
specific angular momentum fluctuations. These accretion 
episodes might lead to the formation of i ntermittent accre¬ 
tion d isks or accretion belts that launch jets (ISchreier & Sokerl 
l2016h . The present study is another step in developing the jit¬ 
tering jets model for CCSNe. This is a preliminary study, and 
the presented results will have to be scrutinized with higher- 
resolution simulations. Also, in future works we will need 
to consider the combined effect of the different sources of 
stochastic angular momentum, as well as the angular momen¬ 
tum contribution from rotation. Magnetic fields may also be 
relevant for the pre-collapse evolution of stochastic angular 
momentum, and they are most definitely important for the 
post-collapse possibility of jet formation. 
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APPENDIX 

A. LOWER RESOLUTION SIMULATIONS 

To check the effects of the simulation resolution on the angular momentum fluctuations, we performed additional simulations 
with the best resolution of the AMR grid (termed flnest-resolution) lower than our nominal simulation. Our computational 
resources prohibit us from increasing the resolution in a full simulation needed when looking at angular momentum of entire 
shells. Other properties of the flow are studied in higher resolution octant simulations. This is presented in AppendixiBl 
Eigure[8]shows the magnitude of the specific angular momentum in different shells for simulations of varying flnest-resolutions. 
The behavior and magnitude are similar, although hard to quantify, due to the stochastic nature of the angular momentum in the 
shells. Eigure|9]show the resolution dependency of the time-averaged value of the specific angular momentum in different shells, 
and of the time-averaged maximal Mach number. As in Eigure|8] the effect of resolution on the specific angular momentum is not 
entirely clear. The maximal Mach number, however, varies very regularly with the changing resolution. Extrapolating linearly to 
infinite resolution reduces the Mach number by a factor of « 2 compared to our nominal simulation. We might expect the angular 
momentum to change accordingly, but this must be studied in simulations of higher resolution, beyond our current computational 
resources. 


B. OCTANT SIMULATIONS 

To further investigate the issue of convergence, we performed several simulations of one eighth of the full simulation with 
reflecting boundary conditions (‘octant’ simulations), with different maximal simulation resolutions. Two simulations, with 
finest resolutions of 800 km and 1250 km, are presented at f = 20000 s in Eigure[T0] Erom Eigure[T0]we can see that the outer 
part of the convective region is similar in the two simulations, while in the inner part the simulation with lower flnest-resolution 
the flow velocity is somewhat higher. We also see by comparing the upper to lower panel on each side that the magnitude of 
the radial velocity component and the magnitude of the circumferential velocity component are similar. This shows that the 
convective flow is well-developed. 

Eigure[TT] shows the kinetic energy due to radial component of the velocity only and due to the circumferential component of 
the velocity only in the convective region, for octant simulations with three different flnest-resolutions, as well as our nominal 
full simulation, at t — 20000 s. The kinetic energy in the outer part seems converged, while the inner part is sensitive to the 
simulation resolution. The radial and circumferential flow components of the kinetic energy are nearly equal, correspondingly 
with the velocity flows seen in EigurefTOl 

Eigure [12] shows the kinetic energy in the region between r = 40000 km and r — 280000 km as a function of simulation time. 
The kinetic energy in the simulations has a tra nsient phase with a large value, and then goes down and plateaus. This behavior is 
similar to that obtained bv IViallet et al.l (1201 3h in their simulations of a convective envelope of a red giant star (see their Eigure 2). 
It can be seen that the simulations are not converged, as the kinetic energy plateau is lower for more highly-resolved simulations, 
as is the time-averaged maximal Mach number. Extrapolating from the values of the maximal Mach number to infinite resolution 
reduces the Mach number by a factor of 2-3 compared to our nominal simulation, similar to the results of the full simulations 
with lower resolution discussed in AppendixiAl 
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Fig. 10.— Presentation of two variables in two simulations in four different quarters. The two simulations are performed in one octant of the full space due to 
limited computer resources. All velocities are in units of cm (so that the highest velocity magnitude presented is 2 X 10® cm s“^). Top: Circumferential 
velocity (i.e., perpendicular to the radial direction) with finest-resolution of 800 km (left) and 1250 km (right). Light shades of orange represent lower velocities, 
and darker shades represent high circumferential flow velocities. Bottom: Radial velocity with finest-resolution of 800 km (left) and 1250 km (right). Red shades 
represent outward motion, and blue shades represent inward motion. The inner oxygen core is omitted from the figure. 
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Fig. 11.— Left: Kinetic energy of the radial velocity component for six different shel ls in octant simulations of varying finest resolution. The values for finest 
resolutions of 800 km and 1250 km are obtained for the simulations shown in Figure [T^ The kinetic energy is multiplied by eight for the presentation of the 
octant simulations. Right: Kinetic energy of the circumferential velocity component (i.e., perpendicular to the radial direction) for six different shells in octant 
simulations of varying finest resolution. The r coordinate is taken in the middle of each shell. The masses of the shells from the inner shell to the outer are 0.25, 
0.25, 0.23, 0.22, 0.19 and 0.17M©. 
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Fig. 12.— Left: Kinetic energy of the convective flow as function of time for our nominal simulation, and for three octant simulations with different resolutions 
(the kinetic energy is multiplied by eight for the presentation of the octant simulations). Right: Average value of the maximal Mach number for t > 3000 s as a 
function of finest grid resolution. The presented linear fit shows the extrapolated value of the Mach number, fA 0.046. 
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